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A  program  of  research  Is  reviewed  that  addressed  transport 
in  physically  small  (submicron-to-ultrasubmicron)  and/or  high 
frequency  semiconductor  devices.  In  devices  of  the  0.25  um 
characteristic  length,  significant  effects  on  the  device  performance 
arise  from  the  fact  that  the  carriers  have  characteristic  response 
times  on  the  order  of  the  device  transit  time.  Effects  due  to 
temporal  (and  spatial)  memory  retardation,  finite-non-zero  collison 
duration,  intra- collisions!  field  effects,  and  quantization  within 
the  device  govern  the  transport.  Additionally,  device  replication 
can  lead  to  synergistic  system  performance.  This  program  of  study 
has  investigated  the  role  these  effects  play  through  utilization  of 
full  two-dimensional  device  simulations  that  incorporate  the  various 
effects. 
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I.  Introduction 


Over  the  last  two  decades,  the  growth  of  large-scale  integration 
(LSI)  has  been  phenomenal,  and  the  application  of  special  integrated 
circuits  and  microprocessors  has  blossomed.  In  fact,  the  complexity, 
as  defined  by  the  number  of  devices  on  an  Individual  integrated  circuit 
chip  has  approximately  doubled  each  year  over  this  time  span  [1].  There 
are,  of  course,  several  factors  leading  to  this  increase  in  complexity, 
including  major  effects  arising  from  increased  die  size,  increased  circuit 
cleverness,  and  reduced  device  size.  There  are  indications  in  recent 
times  that  this  growth  is  slowing  [2],  and  in  fact  Moore  [1]  suggested 
that  this  might  arise  as  a  consequence  of  a  reduction  in  further  increases 
in  circuit  design  cleverness.  It  is  not  at  all  clear  that  this  is  the 
case,  and  it  seems  equally  likely  that  the  reduction  in  growth  is  due  more 
to  the  fact  that  we  are  approaching  some  limitations  on  our  fabrication 
technology  and  on  our  understanding  of  the  physical  limits  of  device 
performance  in  Integrated  systems.  Certainly,  the  current  trend  of  extra¬ 
polation  of  technology  and  design  techniques  can  be  expected  to  be  subject 
to  severe  limitations  as  we  enter  the  very-large-scale  Integration  (VLSI) 
region,  even  though  current  design  technology  has  certainly  achieved  relatively 
small  devices  in  MOS  [3-5],  bipolar  [6],  and  compound  semiconductor  [7] 
technology. 

From  the  above  discussion,  it  is  apparent  that  one  can  reasonably  ask 
just  how  far  can  we  expect  device  size  to  be  reduced,  and  do  we  understand 
the  physical  principles  that  will  govern  device  behavior  in  that  region  of 
smallness.  Previously,  Hoeneisen  and  Mead  [8,9]  discussed  minimal  device 
size  for  MOS  and  bipolar  devices  based  upon  extrapolation  and  scaling  of 
our  present  device  structures.  They  considered  a  general  study  of  the 
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static  limitations  imposed  by  oxide  breakdown,  source-drain  punchthrough, 
ionization  in  the  channel,  etc.,  on  typical  MOS  structures,  and  similar 
effects  in  the  bipolars.  Although  they  suggested  a  minimum  channel  length 

O 

of  1100A  for  the  MOS,  for  example,  it  was  found  that  a  more  typical  device 

o 

would  have  a  channel  length  of  2500A.  Others  have  also  approached  this 
topic  In  a  similar  manner  and  have  also  considered  the  role  of  doping 
fluctuations,  power  dissipation,  noise,  and  interconnections  in  limiting 
packing  density  of  these  devices  [10].  And  yet,  for  all  of  these  studies, 
it  is  evident  that  the  emphasis  has  been  on  other  questions  than  whether  or  not 
our  current  understanding  of  device  physics  is  pertinent  to  devices  on  that 
spatial  scale  [11].  Present  day  electron-beam  lithography  and  photoresist 
technology  is  capable  of  addressing  devices  of  that  size,  and  In  fact  is 

o 

apparently  capable  of  producing  llnewldths  (and  device  sizes)  of  250A  [12]. 

O 

With  modifications  of  resist  structures,  even  80 A  lines  have  been  produced 
[13]. 

o 

Even  though  Hoeneisen  and  Mead  suggested  that  a  size  of  1100-2500A 
channel  length  might  be  the  minimum,  current  research  endeavors  have  pushed 
to  and  even  beyond  these  limits.  Rockwell  [14],  has  made  MOSFETs  of  0.25  pm 
channel  length  and  width  and  operated  them  in  a  ring  counter  configuration. 

Texas  Instruments  has  used  selective  etching  to  make  MOSFETs  of  0.1  ym  channel 
length  [15],  and  IBM  has  fabricated  devices  of  a  few  hundred  Angstrom  channel 
length  [16].  Similar  progress  is  underway  for  MESFETs  as  well,  where  the 
critical  length  remains  the  gate  length.  Granted  that  technological  momentum 
is  pushing  to  ever  smaller  devices,  and  the  technology  is  there  to  prepare 
really  very  small  devices,  it  becomes  obvious  that  we  must  now  ask  whether  our 
physical  understanding  of  devices  and  their  operation  can  indeed  by  extrapolated 


2 


down  to  very  small  space  (and  time)  scales,  or  do  the  underlying  quantum 
electronic  principles  prevent  a  down-scaling  of  the  essential  semi-classical 
concepts  upon  which  our  current  understanding  is  based?  The  purpose  of  the 
current  work  and  the  work  proposed  here  is  to  address  just  this  question. 

In  the  following  sections,  we  review  the  progress  made  under 
current  ARO  grant  DRXR0-15435-EL,  although  part  of  this  is  deferred  until 
Section  III.  The  ARO  program  is  primarily  concerned  with  the  physics  of 
submicron  devices. 


II.  Physics  of  Submicron  Devices 

In  our  research,  we  have  attempted  to  lay  a  conceptual  framework  for  an 
ultimate  physics  of  small  devices  and  the  modeling  necessary  to  characterize 
these  devices.  We  borrow  heavily  from  quantum  theories  of  nonlinear 
electronic  processes.  Some  of  these  have  already  been  illustrated  to  be  of 
importance,  in  hot-electron  transport  for  example,  but  we  recognize  that 
the  totality  of  such  studies  is  small  and  not  fully  developed.  Many  of  the 
effects  will  be  far  more  significant  in  the  lighter  mass  compound  semiconductors 
[17]  than  in  the  silicon  currently  used.  To  approach  this  problem,  we  find  it 
convenient  to  define  two  dimensional  scales  for  small  devices  (where  the 
dimension  might  be  channel  length  for  example):  the  medium-small  device 

O 

(MSD)  with  typical  spatial  extent  of  2500A,  and  the  very- small  device  (VSD) 

O 

with  typical  spatial  extent  of  250A.  The  reader  will  note  that  these  dimensions 
correlate  with  the  typical  small  MOS  device  discussed  by  Hoeneisen  and  Mead 
[8]  and  with  the  small  structures  prepared  by  Broers,  et  al . ,  [12].  MSD's 
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provide  a  scale  upon  which  It  Is  not  unreasonable  to  expect  that  the  physics 
will  be  well  approximated  by  extrapolations  from  present  day  device  physics 
but  with  certain  new  features  anticipated  from  quantum  physics  which  modify 
considerably  the  device  performance.  MSD's  might  therefore  be  similar  in 
concept,  but  not  necessarily  in  performance,  to  present  day  structures. 
However,  on  the  VSD  scale,  a  fully  quantum  mechanical  understanding  appears 
necessary,  and  indeed  new  device  structures  may  well  be  possible  which 
fully  exploit  the  possibilities  provided  by  quantum  electronic  phenomena. 

Here  the  transport  may  well  be  more  ballistic  in  nature  rather  than  statisti¬ 
cal  as  described  by  current  transport  approaches.  This  scale  has  certain 
very  real  advantages  over  present  day  device  scales,  since  some  limiting 
processes  such  as  dielectric  breakdown  (impact  ionization)  will  be  largely 
eliminated,  although  interband  tunneling  will  probably  be  present  but  at 
much  higher  field  strengths. 

In  approaching  the  question,  we  have  introduced  the  archetypal  MSD  as  a 

O 

MOS  device  of  channel  length  2500A  and  discussed  the  relevant  spatial  and 
temporal  scales  that  are  important  for  transport  within  this  structure  [18] 

We  have  examined  the  important  time,  distance,  velocity,  and  energy  scales 
that  are  important  in  small  semiconductor  devices.  For  this  medium-small 
device,  we  find  that  we  are  already  approaching  the  time  and  distance  scales 
that  raise  serious  objections  to  usage  of  a  semi -classical  approach  to 
transport  theory  within  the  device.  As  a  consequence,  major  modifications 
are  required  to  both  the  collision  interactions  in  the  Boltzmann  transport 
equation,  due  to  intra-collisional  field  effects  arising  from  a  breakdown 
of  the  independent  perturbation  approximation,  and  to  the  current  response 
equations  within  the  device.  In  this  last  case,  effects  arise  not  only  due 
to  spatial  and  temporal  relaxation  of  the  carriers  and  charge  density,  but 


also  due  to  retardation  effects  within  the  current  response.  The  MSD 
thus  represents  not  only  the  logical  lower  limit  expected  from  scaling 
arguments,  but  also  the  lower  limit  for  which  performance  can  be  expected 
to  follow  extrapolations  from  current  experience.  Even  here,  though,  we 
expect  major  modifications  in  the  detailed  response  functions  for  the 
device.  As  a  result,  frequency  and  power  limitations  in  small  devices  can 
be  expected  to  be  significantly  different  from  predicted  values  based 
upon  extrapolations  from  large  devices.  Let  us  now  consider  the  VSD. 

Although  it  has  been  demonstrated  that  lithography  is  capable  of 

O 

fabricating  linewidths  and  spacings  of  250A,  this  is  not  a  sufficient 
condition  for  actually  fabricating  a  device.  However,  we  have. demonstrated 
how  such  a  device,  a  p-channel  MESFET,  can  in  fact  actually  be  fabricated 
from  current  technology,  utilizing  two  side  processing  of  the  semiconductor 
wafer,  here  assumed  to  be  silicon  [19].  However,  it  was  pointed  out  that 
already  in  the  MSD,  the  limits  of  applicability  of  the  Boltzmann  transport 
equation  (BTE)  were  reached.  This  implies  that  for  the  VSD,  transport  must 
be  considered  from  a  fully  quantum  mechanical  viewpoint.  The  foundation 
for  this  was  laid  by  establishing  that  a  distribution  function  picture  is 
still  valid.  We  have  carried  this  forward  utilizing  the  full  density  matrix 
for  the  carriers,  phonons,  and  environment.  We  showed  that  a  BTE-like 
equation  results,  although  the  specific  meaning  of  the  effective  distribution 
function  is  obviously  different.  The  role  of  the  environment  is  especially 
important,  since  in  many  cases  it  can  provide  a  renormalization  of  the  energy 
scales  such  as  in  size-quantization  of  superlattices.  In  the  MSD,  it  was 
shown  that  a  displaced  Maxwellian  distribution  function  was  appropriate  [18,20], 
For  the  VSD,  we  have  shown  that  since  a  BTE-like  transport  equation  is 
obtained,  moment  equations  can  be  developed,  as  in  the  MSD,  using  an  analgous 


formulation  to  the  displaced  Maxwellian,  i.e.-a  "displaced"  density  matrix 
form  of  the  distribution  function. 

If  one  can  summarize  the  difficulties  that  arise  in  these  small  devices, 

the  primary  problem  is  the  very  fast  time  scales  inherent  in  small  devices. 

See,  for  example.  Fig.  1,  where  we  plot  the  appropriate  scales  that  one 

can  encounter  in  small  devices.  An  electron  travelling  at  10^  cm/sec  can  cross 

-12 

0.1  micron  length  in  10  sec.  On  this  time  scale,  the  electrons  encountering 
a  high-field  region  of  this  dimension,  such  as  the  pinch-off  region  in  a 
MOSFET,  do  not  have  adequate  time  to  establish  any  sort  of  equilibrium 
distribution,  a  point  made  rather  pointedly  earlier  by  Ruch  [21],  and  by 
Maloney  and  Frey  [22].  Additional  complications  arise  from  the  fact  that 
the  collision  duration  is  no  longer  negligible  on  this  time  scale  and  strongly 
affects  the  transport  dynamics  [23,24,25], 

In  the  following  sections,  we  examine  the  consequences  of  changes  that 
the  short-time  scales  play  on  the  transport  of  carriers  in  these  small  devices. 
First,  we  look  at  the  effect  of  the  non-negl igible  collision  duration  and  the 
concommitant  intra-col li si onal  field  effects.  Then  we  turn  to  its  effect 
on  the  transport  equations,  especially  with  regard  to  transient  dynamic 
response.  In  this  latter  regard,  the  concept  of  diffusion  on  this  time 
scale  must  be  carefully  examined. 

A.  The  Intra-Collisional  Field  Effect  (ICFE) 

The  Boltzmann  transport  equation  (BTE)  has  long  been  the  basis  for 
semi-classical  transport  studies  in  semiconductors  and  other  materials. 

Its  utility  also  stems  from  the  fact  that  is  is  readily  transformable  into 
a  path-variable  form  which  can  be  adapted  to  numerical  solutions  for  compli¬ 
cated  energy-dependent  scattering  processes  [26,27].  In  this  form,  the 
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BTE  is  often  referred  to  as  the  Chambers -Rees  path  integral  equation,  and  serves 
as  the  basis  for  Monte  Carlo  and  iterative  calculation  of  transport.  However, 
the  BTE  is  valid  only  in  the  weak  coupling  limit  under  the  assumptions  that 
the  electric  field  is  weak  and  slowly  varying  at  most,  the  collisions  are 
independent,  and  the  collisions  occur  Instantaneously  in  space  and  time. 

Each  of  these  approximations  can  be  expected  to  be  violated  in  future  sub¬ 
micron  dimensioned  semiconductor  devices.  We  have  shown  that  in  such  devices, 
the  time  scales  are  such  that  collision  durations  are  no  longer  negligible  when 
compared  to  the  relevant  time  scale  upon  which  transport  through  the  device 
occurs  [18,19].  In  this  situation,  even  for  time-independent  fields,  the 
quantum  kinetic  equations  are  nonlocal  in  time  and  momentum.  It  may  be 
recalled  that  the  BTE  can  be  rigorously  derived  from  the  density  matrix 
Liouville  equation  formulation  of  quantum  transport  [24,28].  In  the  BTE 
approach,  the  collision  terms  are  derived  under  the  assumption  that  the  colli¬ 
sions  occur  instantaneously,  which  is  a  reasonable  approximation  when  the 
mean-time  between  collisions  is  large.  At  high  fields,  such  as  will  occur  in 
very  small  devices,  the  collision  duration  is  significant  and  correction 
terms  must  be  generated  (for  the  BTE)  to  account  for  the  non-zero  time 
duration  of  each  collision.  If  the  instantaneous  collision  approximation 
that  leads  to  the  BTE  is  relaxed,  an  additional  field  contribution  appears 
as  a  differential  super-operator  term  (see,  e.g.,  the  discussion  in  ref.  18) 
in  the  collision  integrals  evaluated  in  the  momentum  representation, 
resulting  in  an  intra-collisional  field  effect  [24,25]. 

The  intra-collisional  field  effect  (ICFE)  can  be  partially  understood 
by  referring  to  Fig.  2.  In  the  Boltzmann  case,  the  collision  occurs 
instantaneously,  so  that  the  carrier  enters  the  collision  sphere  at  a  and 
exits  at  b.  However,  the  collision  does  not  occur  instantaneously,  but 
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requires  a  non-zero  collision  duration  tc.  In  this  case,  it  can  now  be 
accelerated  by  the  field  during  the  collision.  Thus  it  exits  not  at  b, 
but  at  b’  some  time  At  =  tc  later.  When  tc  begins  to  become  comparable 
to  t,  the  mean  time  between  collision,  this  ICFE  will  have  a  significant 
effect  on  the  transport  dynamics,  particularly  in  the  transient  response 
region. 

The  mathematical  details  of  the  ICFE  have  previously  been  given  [24], 
so  we  shall  not  go  into  these  details  here.  Rather,  we  shall  give  some 
of  the  supportive  evidence  for  the  observability  of  the  effect.  In  very 
large  fields,  such  as  can  occur  in  Si Og  near  breakdown,  the  ICFE  can 
indeed  be  very  significant  [29],  Two  major  modifications  of  the  scattering 
integral  occur  as  a  result  of  this  intra-col lisional  process.  First,  the 
total  energy-conserving  6-function  is  broadened  by  the  presence  of  the 
electric  field.  Second,  the  threshold  energy  required  for  the  emission 
of  an  optical  phonon  is  modified,  which  causes  a  shift  (in  energy)  of  the 
6 -function  .  This  latter  process  is  easily  understood  in  physical  terms. 

The  argument  of  the  energy  conserving  6-function  is  just 

Ef  -  E1  ±fi«o0  *  E(pf)  -  Etf.)  ±*<o0  (1) 

but  the  initial  and  final  momenta  evolve  during  the  collision  as 
t 

p(t')=p-  |  eE(t")dt"  ,  (2a) 

t'  t 

p'(t')  =  p’  -  eE(t")dt"  .  (2b) 

t' 

In  the  emission  of  an  optical  phonon  ,  where  the  electron  is  scattered 
against  the  electric  field,  the  field  will  absorb  a  portion  of  the  electron's 
energy  during  the  collision,  and  hence  a  reduction  in  energy  loss  to  the 
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lattice  will  be  favored.  The  opposite  effect,  an  enhancement  in  energy 
loss  to  the  lattice,  occurs  for  emission  along  the  electric  field.  These 
effects  can  be  incorporated  into  the  appropriate  scattering  integrals  in 
the  iterative  technique,  and  this  has  been  carried  out  [29].  For  electric 
fields  above  (5-6)  x  10^  V/cm,  the  broadening  and  shift  of  the  scattering 
resonances  produce  a  noticeable  effect  upon  the  velocity-field  relationship 
and  a  reduction  in  threshold  can  be  observed  in  the  calculated  impact 
ionization  rates  in  SiO^.  In  fig.  3,  we  show  these  rates  as  measured  by 
Solomon  and  Klein  [30]  and  as  calculated  earlier  [29].  It  is  necessary 
to  include  effects  arising  from  the  ICFE  to  adequately  fit  theory  to 
experiment.  The  second  major  effect  of  the  ICFE  arises  for  scattering 
between  non-equivalent  sets  of  valleys  [31].  There,  the  ICFE  causes  a 
virtual  lowering  of  the  final  state  valleys  with  respect  to  the  initial  state 
valleys.  However,  this  effect  is  not  large  in  realistic  devices. 

Although  the  ICFE  is  exceedingly  large  in  SiO^  because  of  the  polar 
nature  of  the  phonons  here,  it  is  also  significant  in  the  case  of  Si.  In 
fig.  4,  we  plot  the  collision  duration  versus  field  and  compare  it  with  the 
momentum  relaxation  time  t  *  m*pe,  where  v  -  v^/F  is  the  mobility.  The  effect 
this  has  on  the  velocity- field  curve  is  shown  in  Fig.  5.  Also  shown  is  the 
countering  effect  of  collision  retardation  discussed  in  the  next  section. 

The  ICFE  is  especially  noticeable  at  high  fields,  where  it  essentially 
eliminates  the  scattering  by  the  low  energy  intervalley  phonon  [32].  This 
phonon  Is  already  weakly  coupled  since  scatters  through  a  first-order  inter¬ 
action  [33]  but  is  normally  an  effective  scatterer  at  high  fields. 

B.  The  Retarded  Transport  Equations 

When  the  ICFE  is  included  as  a  modification  of  the  lowest  order  kinetic 
equations,  a  high  field  quantum  kinetic  equation,  which  replaces  the  BTE, 
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Fig.  3.  Ionization  rates  in  Si02  as  a  function  of  field.  The  dashed 

curve  is  for  m,/m  *  100  and  Ignores  the  subthreshold  phonon 
n  e 

emission.  The  solid  curve  takes  into  account  subthreshold 
emission,  for  -  2.9mQ.  Data  is  due  to  Solomon  and  Klein 
[17].  This  is  fairly  substantive  evidence  for  the  presence 
of  the  ICFE. 


Drift  Velocity  (cm/sec.) 


Fig.  5.  The  velocity-field  curve  for  electrons  in  SI  as  calculated 
by  a  drifted  Maxwellian  approach.  At  high  fields,  the 
field  weakens  the  collisions  (dot-dashed  curve,  b)  causing 
an  Increase  of  velocity  although  the  colllslonal  retardation 
stops  this  effect  (solid  curve).  The  low-frequency  curve 
normally  encountered  Is  the  dashed  curve,  a. 
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is  found  as  [24,34] 


_  ^  t 

+  eE(t)  •  Vp  f(p,t)  «  dt‘ |{S(p,p';t,t' )  f(p',t') 

o  p 

■S(p' ,p;t,t* )  f(p,f )}  ,  (3) 

where  the  momenta  p,p'  are  explicit  functions  of  the  retarded  time  t'  on 
the  right-hand-side  through  the  relationship  of  Eq.  (2),  and  the  transition 
terms  S  take  the  form,  for  inelastic  phonon  scattering, 

S(p,p'  ;t,t' )  =  Re  £  (^L)  exp(-  — ~  )  x  (4) 

q  r  t' 

x  (Nq  +  j  +  |V(q)|2  6+^,  +  ^  exp[-i  j  ^  e(p,p';tM)]  , 

t 

where  e  is  the  argument  of  the  normal  6-function  and  is  given  by  Eq.  (1). 

The  two  exponential  factors  in  (4)  are  related  to  the  joint  spectral  density 
function,  which  reduces  to  an  energy  conserving  6-function  in  the  instantaneous 
collision,  low- field  limit,  e(p)  is  the  quasi -particle  renormalized  energy, 

■ti/T r  is  the  joint  linewidth  due  to  collisional  broadening  of  the  initial  and 
final  states,  and  n  takes  the  values  +1,  -1  for  phonon  emission  or  absorption, 
respectively,  in  the  in-scattering  term.  For  the  out-scattering  term,  the 
role  of  p,  p'  are  interchanged,  although  this  does  net  upset  detailed  balance 
in  the  equilibrium  sense. 

In  small  semiconductor  devices,  where  the  dimensional  scale  is  of  the  order 
of  1.0  yin  or  less,  the  carrier  concentration  will  in  general  be  relatively 
high  [8].  Frohlich  [35]  first  pointed  out  that  the  Isotropic  part  of  the 
carrier  distribution  function  is  Maxwellian  provided  that  the  carrier  con¬ 
centration  exceeds  a  certain  critical  concentration;  i .e. --provided  that  the 
rate  of  inter-carrier  energy  exchange  was  sufficiently  large.  Under  these 
conditions,  the  anisotropic  terms  in  the  distribution  function  are  small  and 
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a  parameterized  distribution  function,  often  called  the  displaced  Maxwellian, 
can  be  utilized  [36].  With  this  assumption,  a  hierarchy  of  moment  equations 
can  be  generated  from  the  BTE,  from  which  the  various  parameters  can  be 
determined  (see,  e.g.,  the  discussion  in  [23,24]).  In  this  work,  we  extend 
these  moment  equations  to  include  first-order  effects  arising  from  the 
non-zero  time  duration  of  the  collision.  Under  the  conditions  prescribed 
above  for  the  drifted  Maxwellian  distribution  function,  the  most  general 
nondegenerate  distribution  which  will  be  stationary  under  the  presence 
of  collisions  is  of  the  form  [37] 

f(E)  =  exp[-Bg(E  -  vd  •  p)]  ,  (5) 

where  8„  and  v.  are  functions  of  time  and  B„  *  l/k„T  .  However,  this  form 
6  0  6  b  6 

is  instantaneous  in  time,  and  we  need  rather  a  proper  retarded  function. 

Such  a  function  Is  suggested  by  the  expansion  of  (5)  as 

f ( E )  =  exp(-eeE)[l  +  vd  •  p6e]  (5a) 

and  a  recognition  of  the  form  required  for  the  memory  functional  [34].  To 

achieve  this  form,  we  specify  f(E)  by  the  ansatz  [this  form  is  not  unique, 

but  it  is  not  felt  at  present  that  the  results  are  critically  dependent  on 

the  form  assumed]: 

t 

f (E)  -  J  dt[6(t  -  t)  -vd  (t  -  t)  '  (p  ^)]  fQ  (E,t)  ,  (5b) 

0 

where  the  two  derivatives  are  required,  respectively,  to  preserve  the  dimen¬ 
sionality  and  to  bring  Bfi  into  the  second  expression  in  the  bracket.  In 
the  limit  of  long  times  (5b)  reduces  to  (5a).  Here,  f  *  exp(-BeE). 

Starting  from  the  density-matrix  developed  form  of  the  BTE,  we  have 
shown  previously  that  collision  terms  derived  in  the  normal  case,  but  modified 
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for  intra-colllsional  field  effects  must  be  convolved  with  a  decay  effect 
over  an  effective  collision  duration  [231  Thus,  the  balance  equations  are 
modified  in  a  straightforward  fashion,  although  the  details  are  much  more 
complicated.  This  latter  follows  from  the  role  of  the  intra-collisional  field 
effects,  which  both  broaden  and  shift  the  resonances,  effectively  lengthening 
the  effective  collision  duration  and  weakening  the  effect  of  the  collision 
itself.  If  (3)  is  Laplace  Transformed,  the  moment  equations  can  be  developed 
by  multiplying  by  an  arbitrary  function  integrating  over  the  momentum  with 

(5b),  so  that  the  moment  equations  are  developed  in  the  transform-domain,  and 
then  retransforming  [38].  This  yields  the  energy-balance  equation 


TE  <1  V.»*  *  •  «<!<*>  -  i 

C  0 


exp(-t/tc)<rE(t  - 


where  we  have  recognized  that  <$>  *  3kgTfi/2  when  <p  -  E. 


t)>dT  (6) 

In  addition,  the 


momentum  balance  equation  is 

3V.  *  t  f  -t'/t 

m*lF  =eF-~  J  vd(t"T)  ]  <rm(T-Tt)>e  dx'dt 

c  o  o 

If  tc-k)  and  v^  is  slowly  varying,  then  (7)  reduces  to  the  normal 
3V  . 

m*  IF  =  eF  "  m*vfrm  • 


(7) 


(8) 


and  rm  is  recognized  as  the  rate  of  momentum  relaxation  (  =  1 / xm ) -  The 
various  convolutions  introduce  a  significant  temporal  retardation  in  the 
rate  of  momentum  relaxation. 


C.  Application  to  Silicon 


Calculations  have  been  carried  out  to  ascertain  the  effect  of  the 
colllsional  retardation  in  a  number  of  semiconductors  [20,38].  Here  we 
treat  only  silicon.  The  parameters  used  for  Si  were  the  normal  values: 


m  =  0.91m  and  m.  =  0.19m  .  Intervalley  scattering  was  treated  by  two 

%  q  t  0 

effective  phonons  [32],  a  630  K  zero-order  coupled  phonon  and  a  190  K 

Q 

first-order  coupled  [33]  phonon.  Coupling  constants  of  9x10  eV/cm  and 
5.6  eV,  respectively,  were  taken  for  these  two  modes,  as  they  give  a  good 
description  not  only  of  the  low  field  mobility  as  a  function  of  temperature, 
but  also  to  the  high-field  velocity  in  quasi -two-dimensional  inversion 
layers  [39].  These  values  differ  somewhat  from  recent  iterative  [40] 
and  Monte  Carlo  [41]  calculations,  but  these  latter  treatments  should  be 
interpreted  carefully  as  they  use  an  Improper  treatment  of  the  low-energy 
phonons . 

In  fig.  6  is  shown  the  transient  reponse  of  the  electrons  to  a 
steady  homogeneous  field  of  20  kV/cm.  Although  the  tc  retardation 
acts  to  speed  up  the  response,  the  dominant  factor  is  the  memory  effect 
of  rm  and  vd.  The  collisional  retardation  speeds  up  the  process,  however, 
primarily  due  to  the  effect  of  slowing  down  changes  in  the  effective 
energy  via  collisional  relaxation. 

D.  The  Semiconductor  Device  Model 

In  order  to  evaluate  the  role  of  the  above  transport  effects  on 
actual  device  performance,  we  have  developed  a  two-dimensional  simulation 
of  a  MESFET.  We  have  used  a  submicron  structure,  with  0.25  urn  source-drain 
spacing,  to  put  us  under  the  scaling  limit.  In  this  manner,  the  device 
performance  is  dominated  by  the  carrier  transport  dynamics.  An  exact 
two-dimensional  numerical  analysis  of  a  submicron  n-InAs  Cold-FET  was 
performed  [42].  The  model  used  a  finite  difference  approach  and  an 
algorithm  developed  to  solve  the  coupled  set  of  differential  equations. 

The  program  solved  Poisson's  equation  using  a  FFT  technique  [43,44]. 
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The  electric  field  and  current  density  were  updated  from  the  solution  of 
Poisson's  equation  and  the  continuity  equation  then  solved  to  update  the 
charge  concentration  using  a  fourth-order  predictor-corrector  integration 
routine  [45].  The  currents  are  obtained  at  the  drain,  gate  and  source, 
and  a  test  on  the  overall  consistency  is  made;  a  negative  outcome  starts 
the  process  over  again,  while  a  positive  outcome  accepts  the  results. 

An  n-channel  device  with  parallel  source  and  drain  contacts  to 
eliminate  geometrical  effects  due  to  co-planar  placement  was  used. 

The  bulk  of  the  FET  consisted  of  an  n-InAs  epitaxial  channel  layer  on 

top  of  an  n'substrate.  All  calculations  were  performed  assuming  a 

17  -3  13  -3 

channel  doping  of  5x10  cm  and  a  substrate  doping  of  1x10  cm  . 

The  drift  current  in  the  channel  was  obtained  using  a  field  dependent 

mobility  curve  obtained  previously  for  InAs  at  77  K.  [45].  The  device 

channel  dimension  was  0.25  and  the  active  depth  0.04  ym.  The  gate 

length  was  0.065  pm.  All  calculations  were  normalized  to  the  width 

of  the  device,  taken  to  be  1.0  ym.  Neumann  boundary  conditions  were 

applied  to  all  insulating  surfaces  in  the  device  and  Dirichlet  boundary 

conditions  to  all  the  contacts.  Ohmic  contacts  were  simulated  on  the 

source  and  drain  by  forcing  charge  neutrality. 

InAs  does  not  readily  form  Schottky  barriers  with  many  materials  [47]. 
This  absence  of  a  barrier  in  n-type  material,  with  even  highly  electro¬ 
negative  metals  such  as  gold,  has  been  confirmed  by  measurements  in  the 
ternary  system  InGaAs  [48].  Thus  the  metal  used  and  the  doping  concen¬ 
tration  under  the  gate  must  be  carefully  considered.  Schottky  barriers 
on  covalent  compound  semiconductors  have  been  suggested  as  being  a  function 
of  the  anion  electronegativity  [49],  and  this  suggests  a  candidate  metalli¬ 
zation.  A  material  that  may  show  some  merit  in  producing  a  Schottky 
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barrier  in  InAs  is  ($N)X  [50].  ( SN) x  has  an  electronegativity  about  0.1  eV 

higher  than  gold,  and  hence  should  be  able  to  produce  a  Schottky  barrier 
with  a  barrier  height  of  about  0.025  eV,  adequate  for  the  device  suggested 
here.  Even  lacking  such  a  Schottky  barrier,  a  barrier  can  be  induced  by 
a  blocking  contact  produced  by  a  doping  notch  [51],  Because  of  the 
uncertainty  in  the  actual  barrier  height,  the  calculated  characteristics 
are  referenced  to  an  effective  gate  bias  V'^  =  ^  -  V  . 

In  :fig.  7,  we  show  the  I-V  curves  obtained  from  the  2-D  numerical 
simulation  of  the  subciicron  InAs  Cold-FET.  Although  a  highly  saturated 
drain  current  is  not  obtained,  good  FET  behavior  is  illustrated.  In 
fig.  8,  we  show  the  charge  density  within  the  device. 

The  transconductance  gm  can  be  extracted  from  the  I-V  curves,  and 

the  value  of  gm  obtained  at  Vd  =  0.1  V  is  approximately  800  millimhos. 

The  maximum  gate  capacitance  is  obtained  by  assuming  a  fully  depleted 

condition.  The  value  of  Cg  was  estimated  to  be  2xl0"16f.  From  these 

values  of  gm  and  C^,  an  RC  switching  time-constant  can  be  evaluated  as 

-14 

Cg/gm  which  is  2.5x10  sec.  However,  the  transient  response  is  routinely 
calculated  in  the  program,  and  a  switching  time  of  7.8xl0~14sec  was  obtained. 


This  value  is  longer  than  that  obtained  from  q  and  C 

m  g 


above,  and  can  be 


explained  if  the  gate-to-drain  parasitic  capacitance  dominates  the  switching 


characteristics.  These  values  lead  to  a  power-delay  product  P-t .  on  the  order 
- 1 8 

of  10  joules,  easily  competitive  with  the  best  estimates  for  Josephson- 


junction  logic  gates. 
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